14th Cambridge Workshop on Cool Stars, Stellar Systems, and the Sun 
ASP Conference Series, Vol. ?, 2007 
G. van Belle 



A Splinter Session on the Thorny Problem of Stellar Ages 
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L. N. Jensen^, Patrick A. Young^'^, Andrea Miglio'', Sydney A. Barnes^ 

Abstract. Accurate stellar ages remain one of the most poorly constrained, 
but most desired, astronomical quantities. Here we briefly summarize some 
recent efforts to improve the stellar age scale from a subset of talks from the 
"Stellar Ages" splinter session at the 14th Cambridge Workshop on Cool Stars, 
Stellar Systems, and the Sun. The topics discussed include both the apparent 
successes and alarming discrepancies in using Li depletion to age-date clusters, 
sources of uncertainty in ages due to input physics in evolutionary models, and 
recent results from asteroseismology and gyrochronology. 

1. Introduction 

^^What is the age of [star catalog] [index is probably one of the most 

frustrating astronomical inquiries that one can make, and also one of the most 
quixotic of tasks to tackle. For studies of the evolution of interesting properties 
related to stars (e.g. planeticity, circumstellar disk evolution, etc.), astronomers 
would ideally like a robust means of estimating the age of any star in the sky. 
There are difficulties with age estimation throughout the stellar mass spectrum, 
with some stellar mass and age regimes more amenable to age-dating than others 
due to rapid structural evolution. To make matters worse for the estimation of 
ages for field stars, most useful age diagnostics can manifest a wide range of val- 
ues amongst stars in supposedly coeval clusters. There are age-dating methods 
that are becoming more precise (e.g. Li depletion boundary, gyrochronology, 
asteroseismology), but it remains to be seen how accurate they are when some 
underlying critical physics and parameters (e.g. the solar metal fraction, con- 
vection, etc.) are still matters of active debate. Of course, the "golden spike" 
of the stellar age scale is the Sun, which has its age constrained by lead isotopic 
ages of th e oldest portions of meteorites (the Ca-Al-rich inclusions; 4567.2 ±0.6 
Myr old; lAmelin et alll2002l l. 
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At the Cool Stars I4 meeting held in Pasadena in 2006, a sphnter session 
of talks was held devoted to discussing contemporary issues regarding the ages 
of cool stars. Here we summarize a subsample of those talks (which due to time 
constraints covered only a subsample of age-related topics) and point the reader 
to recent studies relevant to the stellar age scale. The text in the following 
sections reflects the current state of the field in the view of the various authors, 
and no attempt has been made to reconcile their results and opinions. 



2. Stellar Ages: How Far Can We Go? (Barrado) 

A survey in the literature will easily show that, for a given stellar association, 
there is a large range of age estimates. These differences arise from different 
facts. First, we have different sets of theoretical models which might produce 
different answers for the same dataset. In any case, the comparison between 
observational prope rties and those predicted by theory is not an easy task (see 
Stauffer et allll995l ). Moreover, stars are individuals, and they have properties. 



such as starspots, accretion rates, initial conditions, etc., which can alter the 
age estimate. As an example, colors can be modified by surface inhomogeneities 
(IStauffer et al.ll2003l 'l which can affect a cluster's placement in a color-magnitude 



diagram. An additional problem is coevality. It is normally assumed that mem- 
bers of an association are born almost at the same time, but this is not necessar- 
ily true. On the other hand, among the methods we use to estimate ages, both 
primary (upper Main sequence and isochrone fitting, eclipsing binaries, spectral 
features, lithium) and secondary (stellar activity, rotation, group membership), 
there are many pitfalls since there are uncertainties in the input physics. More- 
over, each of these methods is valid for a given range of stellar ages and masses, 
a fact that is not always honored. All these sources of uncertainty make the 
stellar chronology very uncertain although sometimes we can find formal errors 
of tenths of Myr, which are far from our present capabilities. 

From my point of view, we should aim for consistency. A given property 
(or properties) allow us to sort a set of stellar associations from the youngest 
to the oldest, even if we can not derive absolute ages. Even if other properties 
produce a different age scale, at least we should be able to have the same kind 
of sequence, with Taurus being younger than the A Ori cluster, which is younger 
than IC2391, and still much younger than the Pleiades. Eventually, we should 
be able to construct different ages scales which should be consistent with each 
other, and should produce absolute and well as relative ages. This only can be 
achieved with a considerable effort in the theoretical and observational sides of 
the problem. Specifically, we should gather very large datasets with the same 
setup for a large number of stellar association, in order to avoid biases. 



3. The Li Depletion Boundary in the /? Pic Group (Jensen) 

Pre-main-sequence (PMS) Li depletion has the potential t o be a valuab le age 
indicator for late-type PMS stars, (see e.g., the review by lJeffriesll200d ). Like 



any age indicator, however, the models that tie a particular Li abundance and 
effective temperature {T^s) to an age must be calibrated with empirical data. 

We observed twelve late-type members of the /? Pic moving group (BPMG) 
at spectral resolution R = 40, 000 using the echelle spectrograph on the 4-meter 
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D'An:ora & Wazzitelli ;i997/-99S;: Siess et al, (2000) Baraffe et al. (199B) 
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Figure 1. Toff vs. Li abundance for stars in the 10-12 Myr-old /3 Pic moving 
group. For each model, isochrones for 1, 3, 10, 30, and 100 Myr are plotted. 



Blanco telescope at CTIO. Spectral types were determined from the TiO band in - 



dices, and converted to T^g using the intermediate scale giv en bvlLuhman I (119991) 



for th e M stars and the main-sequence (MS) scale from iKenvon &: Hartmann 



( 19951 ) for earlier types. We determined Li abundances by comparing the ob- 



served spectra to a grid of high-resolution synthetic spectra. 

Fig. [J shows Li abundance data for t he twelve stars overlaid on three sets 
of m odels: (jP'Antona k Mazzitellil Il997l . Il998l : ISiess et al]l2000l : iBaraffe et al 



19981 ). These stars all have ages of ~ 10-20 Myr based on their HR diagram 



positions, but the systematic trend observed in their Li depletion is that full 
depletion of Li occurs much more rapidly for the M stars than the models pre- 
dict. In each of the models, the right side of the dip in the lithium depletion 
isochrones, which indicates the lithium depletion boundary (LDB), consistently 
falls at effective temperatures too high for our data if we assume the HR diagram 
ages are correct. Three stars with spectral types of M4 show fully-depleted Li, 
with an upper limit of log A{Li) < —1; none of the models predicts significant 
Li depletion in stars of this Tgfj until ages > 30 Myr. 

Some authors have questioned whether the BPMG is truly a coeval group 
of stars with a common origin. Even if one accepts an age spread in the sam- 
ple, however, the components of individual binary or multiple systems should 
be coeval with each other. Examination of the Li abundances in individual 
systems in the sample shows the same discrepancies as the sample as a whole; 
the pair BD-17 6128 A/B, and the triple systems HD 155555 A/B/C and GJ 
803/799N/799S, are inconsistent with having a single Li-depletion age for all 
components. The direction of the systematic discrepancies reported here are 
consistent with previous work showing that the LDB occurs at a cooler T^s for 
a given age than the models predic t. Our results are also consistent with the 
work of White &: Hillenbrand ( 20051 ). who suggest that the M3 classical T Tauri 



star could be a young member of Tau-Aur despite its large Li depletion. 
4. Lithium Depletion Ages and Age Spreads (Randich) 

HR diagrams of several young clusters and associations suggest the presence of 
age spreads among cluster members. Understanding whether these spreads are 
intrinsic or rather due to errors in the determination of stellar parameters and 
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misplacement of the stars on the HR diagram, is crucial for gaining insight into 
the cluster star formation (SF) histories and, in particular, to put constraints on 
the overall duration of the SF process (Hillenbrand, this volume). Constraints 
on aspects of SF and early (sub)stellar evolution can also be gathered through 
the investigation of PMS clusters with ages ~ 10-50 Myr. SF in PMS clusters 
has ceased, implying that their members have attained their final masses and 
that these clusters do not present the complications that are intrinsic to SF 
environments. On the other hand, PMS clusters are generally too young to be 
dynamically relaxed and thus reveal a cleared view of the final product of the SF 
process. In this context, using the multiplex instrument FLAMES on VLT/UT2, 
we have obtained medium/high resolution spectra of low-mass members of the 
Orion Nebula Cluster (ONC), the a Ori cluster, and the PMS cluster IC 4665 
(|de Wit et al.ll2006l 1. Our main goal was the measurement of their lithium con- 
tent, which provides an independent clock to derive the ages and age spreads. 
ONC and a Ori - (jPalla et al.ll2005l . l2007l : ISacco et aDl2007l ): We have 
observed 141 high-probability cluster members (membership probability from 
proper motion >90%) down to -^O.IMq in the ONC and 98 cluster candidates 
in a Ori. Li abundances were derived from veiling corrected pseudo-equivalent 
widths (pEWs) of the Li i 670.8 nm feature. Radial velocities were also measured, 
along with the strength of the He 667.8 nm (ONC) and Ha (cr Ori) emission 
lines. All the ONC targets were confirmed as members, while in a Ori 59 bona 
fide members were identified. Li abundances are consistent with the interstellar 
value for most of the sample stars in both clusters. However, six stars in the 
ONC and three stars in cr Ori show Li depletion factors between ~2-50 (ONC) 
and >100 (cr Ori). Nuclear ages and masses for four Li-depleted stars in the 
ONC and two of the three Li-depleted stars in a Ori are in excellent agreement 
with isochronal ages and masses. Liferred ages are in the range ~10-30Myr, 
much older than the bulk of the cluster populations. Hence, for both regions, 
there is independent evidence for a population of old stars mixed with the much 
larger assembly of young objects, supporting the suggestion that SF in the two 
clusters has not occurred in a short single burst, but has rather continued for a 
long time, much in excess of the dynamical time scale. 

IC 4665 — (Manzi et al. in prep.): With a turn-off age of 36 Myr, this 
cluster represents a good candidate PMS cluster. We observed 158 candidate 
members and confirmed membership for 45 of them based on radial velocity, 
presence of Ha and/or Li absorption line (depending on mass). The Im vs. Im-z 
color-magnitude diagram of the confirmed members is shown in Fig. O together 
with isoc hrones and the pre dicted location of the LDB from the evolutionary 
tracks of Bar affe et al. ( 19981 ). The figure shows a sharp boundary: the faintest 
star without a detected Li line has Im = 16.39, while the brightest star with Li 
has 1=16.61. The LDB was determined as the mean of the two values, assuming 
d = 370 pc, and reddening Kj^ = 0.326. After considering all of the observational 
uncertainties, this yields Mj^, = 8.52 it 0.40. Accordingly, IC 4665 has an LDB 
age of fLDB = 27 =b 5 Myr. This value confirms the PMS nature of IC 4665 and 
agrees with the turn-off age, the first secure instance for a PMS cluster. 
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Figure 2. Im vs. (Im — diagram for members of IC 4665. Filled and open 
symbols indicate stars with and without Li. The 100 Myr ( solid line) and 30 
Myr (short-dashed line) isochrones from lBaraffe et aT] (|1998i) are also shown, 
along with the predicted LDBs for 20, 30 and 40 Myr. 

5. Theoretical Uncertainties in Evolutionary Tracks and Ages (Young) 

The accuracy of stellar ages determined from fitting to evolutionary tracks or 
isochrones is limited by the accuracy of the tracks themselves. Inaccuracies 
can arise from missing or incorrect physics in the models as well as physical 
quantities known with limited precision, such as metallicity. Some of the most 
important sources of error are discussed here individually, but it is important to 
remember that there are cross terms between all of these quantities. 

The most important problem is the hydrodynamics in stellar models. Be- 
cause evolution calculations are treated as a series of snapshots of static states, 
dynamics are often ignored. Some elements such as rotation have been added, 
but are generally treated analytically and in isolation. Deficiencies are often 
corrected by the inclusion of free parameters calibrated by observations. A cal- 
ibration appropriate for one stellar type is generally incorrect for others. As a 
result, mixing in stars in very poorly described and not predictive. 

Convection is defined by thermodynamic gradients, which is appropriate 
only for the onset of convection. Once convection is established, the dynamically 
unstable part of the star is larger than the thermodynamically unstable region, 
and other processes besides convection can contribute to mixing. Evaluating the 
stability of the stellar stratification against all the available sources of kinetic 
energy, including turbulent convection, waves, and bulk rotational flows gives 
a self-consistent and predictive description of the hydrodynamics. Differences 
from standard models can be large, especially for massive stars. Larger mixed 
regions result in higher luminosities, larger radii, efficient angular momentum 
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transport, and transport of chemical species to the photosphere. On the MS 
model ages for a sample of eclipsing binaries from 1 to 23 Mq change by 15-400%. 
Radiative transport can be important for stars with deep convective envelopes. 
The diffusion time of energy from a fluid parcel can become comparable or even 
less than the crossing time of the convection zone. Convection becomes very 
inefficient, leading to larger radii, and ages can differ by lO's of %. 

Even complete and correctly implemented physics can cause problems due 
to the inherent lack of precis ion in our determination of physical quantities. A 
comparison of the OPAL and[T immes fc ArnettI (|l999l ) equations of state (EOS) 
gives surprising results. The two formulations differ by at most 2% under any of 
the conditions tested, yet they produce very different radii on the Hayashi track. 
Ages determined from each model can vary by ~60% or so. The main sequences 
are not reconcilable without changing the composition or convective treatment. 
This is likely the largest source of uncertainty for sub-solar mass tracks. 

Composition, unsurprisingly, plays a significant role. For a 10% change in 
composition, deeply convective stars can have ages different by ~10%, while 
MS ages can change by factors of several. Given that determinations of the 
solar composition have changed by roughly a third in five years, this is a serious 
problem, especially when precise metallicities are unavailable for most stars. 

The conversion to observables adds another layer of complexity beyond the 
scope of this discussion. Model atmospheres, temperature to color conversions, 
rotational evolution (which require magnetohydro dynamics to model), and cor- 
rect distances are all necessary. 

So far we have discussed these sources in isolation. When matching model 
tracks, radius is the primary comparison, whether directly or through T(.ff. 
Individual effects drive radius in opposite directions, and so cannot be calibrated 
out. As a very simplified example of this interplay, consider a change in the 
equation of state. This changes the stratification of the star. A change in the 
stratification changes the extent and vigor of convection. Hydrodynamics enters 
to change the depth of convection and introduces new mixing processes that 
change the compositional structure and thus stratification of the star. These 
processes also transport angular momentum, which feeds back on the stability 
and mixing. These change shift opacities as well as solutions for the EOS. The 
problem must be solved self-consistently in order to reduce the uncertainties and 
produce an accurate result. 

In short, hydrodynamics introduces changes of lO's to lOO's of % in age 
determinations. Radiation transport adds lO's of % on the Hayashi track. The 
EOS is a surprisingly large contributor at lO's of % at 1 Mq. Composition can 
be a large effect for relatively small changes. All of these sources influence each 
other. These cross terms can result in much larger synergistic errors, or, more in- 
sidiously, the right answers for the wrong reasons. Fortunately, multidimensional 
simulations are allowing us to understand and generalize the hydrodynamic pro- 
cesses in stars, resulting in a considerable improvement in the models. 



6. Asteroseismological Constraints on Ages (Miglio) 



Though the potential of solar-like oscillations as an effecti ve tool to determine 
stellar ages has been known for years (e.g. Gough 19951 ). it is just recently 
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that, thanks to the improvement of observational techniques (see Bedding, this 
volume) , the seismic modeling of bright nearby stars has become reality. 

The uniqueness of a seismic inference of age resides in the fact that oscilla- 
tion frequencies carry direct information on the structure of the central regions 
of the star. A well-known indicator of stellar age for MS stars is the so-called 
small frequency se paration (Jz^) which, as shown by the asymptotic theory of 
stellar oscillations ( Tassoullll980l ). represents a probe of the sound speed in the 
energy-generating core. The evolutionary state of a s tar can also be inferred 
by ot her suitable combinations of frequencies (see e.g. [Roxburgh Sz Vorontsov 
20031 1 or, in the case of models of evolved stars, by the detection in the os- 
cillation spectrum of modes of mixed pressure an d gravity character (see e.g. 
Christensen-Dalsgaard. Bedding, fc Kieldsenlll995l ). 

As an example of one of the first attempts to model solar-like stars in- 
cluding seismic constraints, we consider the visual binary system a Cen that, 
thanks to the combination of precise "classical" constraints (masses, radii, chem- 
ical composit ion, luminosity) and the detection of solar-like osci l lations in both 
components (Bouchv &: Carrieil l2002l : [Carrier BourbanI l2003l : iBedding et al 



2004 : iKieldsen et al-y gOOSil . has recently been the subject of several theoreti 



cal studies (see e.g. Thevenin et all l2002l : iThoul et"al 20031: Eggenberger et al 
20041 : iMiglio fc MontalbanI 120051 '). As predicted bv lBrown et al.l (|l994l '). the in 



elusion of in the modeling of the system allows one to significantly reduce the 
estimated uncertainty on the age of a Cen. If the asteroseismic constraints are 
not included in the fit, a value of 8.9 it 1.9 Gyr is found, whereas the inclusion of 
seismic cons traints red uces the estimated age (and its uncertainty) to 5.8 ± 0.2 
Gyr (see Mig lio fc Mo ntalban 200^). 

The accuracy of such age estimates is limited, however, by the poor fre- 
quency resolution of current seismic data and, more importantly, by uncertain- 
ties in the physics included in the stellar models. Nonetheless, once precise 
seismic data will be available in observationally well-constrained targets (e.g. 
visual binaries), asteroseismology will also be able to go beyond the determi- 
nation of global parameters, and directly test stellar models - not only adding 
precision, but also accuracy to the age estimates. 



7. Gyrochronology (Barnes) 

That the rotation of cool stars slows as a fu nction of age was first put on a firm 
observational footing bv ISkumanich ( 19721 ). The wealth of observational data 
for cool stars in clusters with "known" ages gathered over the past decade can be 
used to construct a useful age estimator from measured stell ar rotat i on pe riods 
{^'gyrochronology"). More details about this can be found in Barnesl ( 20031 ) and 
Barnes (submitted to ApJ). 

Beyond an age of ~100Myr, the rotation rate of any late-F to early-M star 
(those with surface convection zones and interior radiative zones) is describable 
by P = f{B—V)xg(t), where P is the rotation period of the star, and / and g are 
separable functions of stellar color, B—V, and age, t, respectively. Both f{B—V) 
and g{t) can be determined purely from the data. From the open cluster data, 
one can show that f{B-V) = (0.77 3±0.011) x (.5-^ -0.4)°-60i^0-024-), g^t) ^an 



be shown to be roughly Skumanich ( SkumanichI 119721 '). and by requiring that the 



solar rotation rate be reproduced at solar age, one finds that g{t) = t' 



0.519±0.007 
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The preceding considerations allow us to write the "gyro age" of a star as: 
log{tgyro) = \ {logiP) - log{a) -bx log{B -V - 0.4)}, where n = 0.519 ± 0.007, 
a = 0.773±0.011, and h = 0.601ib0.024. One can also derive an error on the gyro 

age. The fractional error is given by f = 2% x A.b + \{lntY + 2P0-6 + (M)2^ 
where x = B — V — 0.4. For late-F to early-M stars of ~1 Gyr, this works out 
to be ~15%. 

This method satisfies all five criteria desired for an age indicator: of finding a 
variable that changes regularly and sensitively with age, of being able to calibrate 
that variable against a very well-defined age, or set thereof, of being able to 
identify and measure the functional form of the variable against its dependent 
variables (this one is even separable), of being able to invert the dependence to 
find the dependence of age on the other variables, and finally, of being able to 
calculate the error as a function of the variables. 
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